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For the Run 2 of the LHC next-to-leading order electroweak corrections will play an 
important role. Even though they are typically moderate at the level of total cross 
sections they can lead to substantial deviations in the shapes of distributions. In 
particular for new physics searches but also for a precise determination of Standard 
Model observables their inclusion in the theoretical predictions is mandatory for a 
reliable estimation of the Standard Model contribution. In this article we review the 
status and recent developments in electroweak calculations and their automation 
for LHC processes. We discuss general issues and properties of NLO electroweak 
corrections and present some examples, including the full calculation of the NLO 
corrections to the production of a W-boson in association with two jets computed 
using GoSam interfaced to MadDipole. 
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I. INTRODUCTION 


The Large Hadron Collider (LHC) was designed to provide a deeper insight into the mi¬ 
croscopic world and increase our knowledge about the fundamental interactions occurring in 
nature. 

The recent discovery of the Higgs boson 1, ij supports the predictions made by the Standard 
Model of particle physics 3-7]. It is not yet clear though, whether the discovered boson is 
exactly the one predicted by the Standard Model, or whether it originates from supersymmetric 
theories, or whether it is something completely different mimicking the role and the properties 
of a Higgs boson. 

There is the hope that the run 2 at the LHC will show not only more phenomenological 
properties of the Higgs particle, but also additional new physics beyond the Standard Model 
(BSM). As in principle it is not known how new physics will manifest itself, one has to look for 
even small deviations from the Standard Model predictions. Completely different new physics 
scenarios can lead to very similar predictions, which makes it difficult to distinguish among the 
models. Besides precise experimental measurements this requires precise theoretical predictions. 

At high energy hadron colliders, final states typically involve hundreds of particles (hadrons, 
leptons, neutrinos and photons). General purpose Monte Carlo (MC) event generators like 


n 


HERWIG |g], HERWIG++ |9|, PYTHIA 6 JLf|, PYTHIA 8 n| and SHERPA ^ can 
provide a fully exclusive description of such complicated reactions. These tools perform several 
tasks ranging from the computation of hard scattering processes to the simulation of multiple 
soft and/or collinear QCD and QED radiation by means of parton shower algorithms. MC 
event generators also provide a description of the non-perturbative effects like hadronization 
and formation of the underlying event. General purpose MC generators play a crucial role 
in collider physics and indeed there is no experimental analysis or realistic phenomenological 
study which does not rely at least in part on these codes. 

One of the drawbacks of general purpose MC event generators is the fact that hard scattering 
processes (and, as a consequence, overall normalizations) are computed at tree-level accuracy, 
while precise theoretical predictions for several processes of phenomenological interest require 
at least the calculation of next-to-leading order (NLO) contributions in perturbation theory. 
Typically, at hadron colliders the largest contributions come from QCD corrections. In the 
last recent years, lots of progress has been made in calculating NLO QCD corrections: new 
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algorithms, such as the ones based on the OPP method 13H15I] or on the generalized unitarity 


techniques 


16 


18j] , have been proposed to compute one loop virtual amplitudes numerically or 


semi-analytically, i.e. bypassing the explicit calculation of one loop Feynman diagrams from 
the Feynman rules of the theory. In particular, a lot of effort has been put into the automa¬ 
tion of such intricate and time consuming calculations and several tools have been developed 


either based on generalized unitarity methods, such as BlackHat [19] and Njet 


on OPP-inspired techniques, such as GoSam 
MadGraph 


25 


261 ]. In codes like OpenLoops 


22], HELAC-NLO 


27J and RECOLA 


20 


231 ]. MadLoop 


or 


24]/ 


one loop amplitudes 


are computed numerically from tree level amplitudes by means of recursion relations applied 
to Feynman diagrams and off-shell currents, respectively. The fact that calculations were not 
carried out any more on a pro cess-by-pro cess basis but obtained by automating the underlying 
building blocks and interfacing them in a standardized way, played a major role in what has 
become known as the NLO revolution. 

The NLO revolution of QCD corrections set a new standard for the data/theory comparison 
at run I of the LHC. While, at the TEVATRON, data were typically compared to the predictions 
of LO multi-parton event generators such as ALPGEN 29], SHERPA and MadGraph that 
were used to generate samples with different hard parton multiplicities merged to parton showers 


in the CKKW [30|, |31| or MLM 32] framework, the standard theoretical benchmark at the run 
I of the LHC for most of the analyses are NLO QCD computations consistently matched with 
parton showers within the POWHEG 33-351 or the MCRNLO 36[_ framework. Methods like 


FxFx 


33], MEPSONLO 38|, M\- MiNLO (40,(411 and UNLOPS (42] were also developed to 


merge parton showers and NLO samples with different parton multiplicity. Both matched and 
merged calculations rely on MC event generators like HERWIG, PYTHIA and SHERPA for 
what concerns showering and non perturbative effects for the description of multiple photon 
radiation off final state leptons the code PHOTOS 43] is also widely used). Needless to say, 
this NLO revolution pushed the frontier of fixed order QCD corrections to the NNLO level of 
accuracy. 

A lot of work has been done in the field of automation of one loop QCD corrections. However, 
as we will point out in the next paragraphs, in view of the 14 TeV run of the LHC there are the 
physical motivations to extend the technical know-how of these automated approaches also to 
the context of electroweak corrections. Because of this point, we decided to focus the present 




































































5 


review on the practical issues that one has to face in order to build a tool for the automated 
computation of NLO electroweak corrections to LHC processes starting from already available 
codes developed for the computation of the different building blocks of an NLO calculation. 


In particular, we will show how this task can be accomplished in GoSam+MadDipole 


44 


46] framework for the computation of the relevant amplitudes and the subtraction of infrared 


singularities respectively, once the effect of electroweak renormalization has been properly taken 
into account. 

The paper is organized as follows. In subsection II Al we give a short overview about the 
phenomenological motivations for the calculation of NLO electroweak corrections at high en¬ 
ergy hadron colliders. In this context, in subsection II Bl we explain the main reasons why the 
automation of 0(a ) corrections will play a crucial role at the 14 TeV run of the LHC and, 
of course, at future colliders. We summarize the state of the art for the automation of NLO 
electroweak corrections, describing very briefly the existing tools and the phenomenological 
studies obtained with such codes. As stated above, with the present review we want to show 
how the automation of NLO electroweak corrections can be obtained starting from already 
existing tools: in particular we work in the GoSam+MadDipole framework. In section UTI we 
discuss some of the issues related to the computation of virtual one loop electroweak correc¬ 
tions, such as the UV and IR regularization scheme dependence, the inclusion of electroweak 
renormalization and the subtraction of infrared singularities. As codes like GoSam work in 
dimensional reduction (DRED), while electroweak renormalization counterterms and infrared 
subtraction terms are usually derived in conventional dimensional regularization (CDR), we 
give the transition rules from CDR to DRED regularization schemes for both renormalization 
counterterms and infrared subtraction terms. We show how these scheme conversion works 
with a few simple examples, checking explicitly the regularization scheme independence of the 
sum of virtual one loop corrections, renormalization counterterms and integrated dipoles. Real 
radiation effects, and in particular the treatment of photon radiation off massless fermions, are 
described in section IIIII As a realistic example of the technical issues described in sections [Hi 
and IIIIl we describe in detail a calculation of the one loop electroweak corrections to the process 
pp —>■ W + jj in the GoSam+MadDipole framework in section IIVI We show the results of 
some internal consistency checks concerning poles cancellation and we briefly discuss the phe¬ 
nomenological impact of the electroweak corrections to the W production in association with 
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2 jets. 


A. Physical motivations 


In the next run of the LHC, with higher center of mass energy and also higher luminosity, 
electroweak radiative corrections will play a more and more important role. 

First of all, the high luminosity will allow to perform electroweak precision measurements. 
This is for example the case for the W boson mass measurements in the charged Drell-Yan 
process, where Mw is determined indirectly from a template fit of the distributions of transverse 
mass Mt = \/2p l T p^(l — cosA^iQ of the lepton neutrino system, of the the lepton px and 
of the missing transverse energy. It is well known from the TEVATRON experience that 
a precise knowledge of the electroweak radiative corrections to the Drell-Yan is mandatory, 
as they alter the shapes of the above-mentioned distributions and thus affect the W mass 
determination. This is the reason why different groups spent a lot of effort to compute the one 


ions to the Drel 


-Yan and several simulation codes have been provic 


55 


ed, 


56] 


loop electroweak correc 

such as HORACE 47, 48], RADY 49-51], SANC 52, 531], WGRAD [54|, WINHAC 
and ZGRAD js7] . All of these codes have LO QCD accuracy and NLO EW accuracy. In 
particular SANC also include NLO QCD corrections, while HORACE includes the effect of 
all order photon radiation properly matched to the 0(a) corrections. NLO EW corrections are 
added to the 0(a 2 s ) ones in the FEWZ code 58], while factorized NLO EW and NLO QCD 
corrections to the single W and Z production matched with QED and QCD parton shower 


have been implemented in the POWHEG-BOX Monte Carlo event generator 


59 


61] 


The higher luminosity will also improve the experimental accuracy for several other processes 
besides the Drell-Yan process. In order to match the foreseen experimental precision, theoretical 
predictions for many processes should include at least the effect of O(a) corrections, as can 
be seen for example by looking at the Les Houches wish-list in Tables 1-3 of Ref. 62]. It 
is interesting to noting how the Les Houches wish-list also points out that, in order to fully 
exploit the physics potential of the LHC, fixed order EW computations should be properly 


matched with parton showers. At present, the issue of matching IN 
parton showers has only been addressed for Drell-Yan in Refs. [59j, 


Higgs decay into four leptons in Ref. 


63] 


LO EW corrections with 
611 ] and very recently for 
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In the context of the high precision determination of the Standard Model parameters such 


as, for example, the W boson mass, the theoretical uncertainties related to parton distri 
functions (PDFs) will eventually become a limitation. Recently new PDF sets 


64 


jution 


671 ] have 


been released which include the constraints coming from LHC run I data. PDFs are obtained 
by fitting experimental data with theoretical computations that are performed at the high¬ 
est available perturbative order which at present is typically NNLO QCD. In Ref. 


electroweak corrections have been considered in the fit of Drell-Yan data 


65J NLO 


68| : the systematic 


inclusion of one loop EW corrections in PDF fits will move in the direction of reducing further 
the theoretical uncertainty related to PDFs. Besides the impact of EW corrections on PDF 
fits, electroweak corrections also affect the description of PDF evolution: QED corrections, for 
instance, can be included in the evolution equations, leading to the definition of the photon 


PDF 


69 


73]. In turn, the knowledge of the photon PDF allows to compute the so-called pho¬ 


ton induced processes (i.e. partonic channels with photons in the initial state) which should be 
consistently included in the computation of O(a) correction s 124211 . At variance with the case 
of QED, the effect of weak corrections to PDF evolution [m], Q has still to be investigated. 

If on one hand the higher luminosity will improve the experimental accuracy, on the other 
hand the increase of the center of mass energy will allow to probe more and more extreme 
regions of the phase space, such as the ones characterized by a large jet multiplicity at high px 
and missing transverse energy (also with additional hard isolated leptons) or by the production 
of boosted objects in association with jets, that are of great interest for both the direct new 
physics searches and also for quite challenging Standard Model measurements (see for instance 


the results of Refs. 


76-81] for W+ multi-jet and Z+ multi 


of direct new physics searches at the LHC we refer to 82-85]). In these extreme kinematical 


et production, as possible examples 


configurations one loop electroweak corrections are dominated by double and single logarithms 
of the typical energy scale of the process over the gauge boson masses (the so called Sudakov 


logarithms 


86 


87]). The Sudakov logarithms basically correspond to the infrared regime of the 


virtual one loop weak corrections, where the gauge boson masses can be neglected if compared 
to the other energy scales involved. In the Sudakov limit the weak boson masses play the 
role of physical regulators for the infrared singularities of the loop diagrams, which in turn 
manifest themselves as logarithms of the energy scale over the infrared cutoffs Mw and Mz 


88h96| . The Sudakov corrections are negative and usually have a moderate impact on the total 































cross sections. However these corrections are sizable in particular in the tails of the pt or Ht 


distributions 


62 


97j where they can be of the order of several tens of percent (—40%, —50%) at 


the 14 TeV LHC and become even larger at the center of mass ener gies of the proposed future 


hadron colliders such as the High-Energy LHC or the hh-FCC 


98 


99j. These large negative 


effects should be included in the theoretical predictions and are a clear indication of the role 
that the O(a) corrections will play at high energy hadron colliders (also for processes with 
many final state particles). 

B. Electroweak corrections for the LHC: quest for automation 


At present full one loop electroweak corrections are available only for a limited cl ass of 


LHC processe s, such a s charged a nd neutra l D rell-Yan. V + 1 


dilepton-l-jets 104j - ll06 ]. single top 107H109 ]. tt 


110 


jet 

V = 

117- 

119 


(including the Z/W decay products), H production in ve ctor boson fusion 
on-shell vec tor bosons or in pole approximation) 


122 


WWZ 


132], H —* 4/ 63|, 133], HA pro duction 134]. tt + H 


124 


129 ]. WW + 1 jet 130 ]. WZZ 


O(a) corrections to Z(—>11) + 2 jets 137], W + n jets (n < 3) 138 ] and tt + H/Z/W 135 . 139 ]. 


135 


r + H 


100- 

1031 

120 

121 


123J, W 


with 


lain. 


136]. More recently, exac t 


Except for the last three examples, all the above computations have been carried out on a 
process-specific approach and are basically limited to processes with four or at most five external 
particles: the main reason of this limitation is the fact that the computational complexity grows 
considerably as the number of external leg particles increases. 

As already stressed in the previous paragraphs, in order to fully exploit the physics potential 
of the second run of the LHC, the computation of the NLO electroweak corrections will be 
mandatory for a very large number of processes (also with large final state multiplicities). It 
should be clear that in this context a process-specific approach for NLO EW calculations is no 
longer feasible and more general and process-independent techniques should be developed. 

As far as only the Sudakov regime is concerned, the universality of the infrared limit of 


weak corrections can be exploited to develop general algorit 


ims for the calculation of the EW 


corrections in the logarithmic 
dakov corrections to diboson 
H 


146 


approximatio n 


93 


95 


96 


124, 


141 


140( | . Following this app roach, th e Su - 


143], vector boson plus multi-jets |97 M144j ], tf+jets [145], 


147] and iL+jet 144 ] production have been computed pointing out further the phe- 
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nomenological impact of the EW corrections at high energies. Order a corrections to dijet, 
Drcll- Yan and it production have been recently included in the MCFM Monte Carlo pro¬ 
gram 148|: both the Sudakov approximation and the full one loop corrections have been 


implemented in order to provide a tool for the fast evaluation of the approximated 0(a) cor¬ 


rections that also allows to asses the validity of the approximated results 


1491. 


Although the Sudakov corrections are universal, this is not the case for the full NLO elec- 
troweak corrections that are supposed to work also for energy scales well below the Sudakov 
regime. Indeed, the non-logarithmic parts of the corrections are strictly process dependent. 
Despite the fact that 0(a) corrections are not universal, what could be cast in a process- 
independent form is the way these corrections are computed, namely the algorithms used to 
compute virtual loop diagrams, to perform the subtraction of IR singularities, to include ultra¬ 
violet renormalization and to compute real corrections. This is the first mandatory step in the 
direction of the automation of NLO EW calculations. 

The development of process-independent algorithms to compute radiative corrections was 
the starting point also for the automation of QCD NLO corrections. Even if in principle the 
strategies developed for one loop QCD calculations can be extended to cover also the case of EW 
corrections, this generalization is not trivial because of several specific aspects of electroweak 
interactions, such as for example the chiral structure of weak interactions (which requires a 
consistent treatment of y 5 matrices in loop diagrams), the treatment of unstable particles (that 
leads to the introduction of complex masses and -as result- of complex couplings), the highly 
non trivial interplay between QED and QCD-like singularities or the presence of very different 
mass scales in loop diagrams. On the real emission side the presence of several mass scales and 
unstable particles typically leads to a phase space with a more involved structure containing 
Breit-Wigner peaks compared for instance to massless QCD. 

Quite recently the issue of automation of one loop electroweak corrections has been ad¬ 
dressed by several collaborations and tools like MadGraph5_aMC@NLO, OpenLoops and 
RECOLA have been developed and applied to perform phenomenological studies on the im¬ 
pact of the NLO EW corrections to the processes pp —>• it + {H, Z, W}, pp W + n jets 
(n = 1, 2, 3) and pp —> Z(U) + 2 jets, respectively. 

In Ref. [138] the factorized NLO QCD+EW corrections to the processes pp —> fU + +l, 2 
and 3 jets have been computed using OpenLoops interfaced with MUNICH and SHERPA. 
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OpenLoops computes LO and NLO amplitudes by means of four-dimensional recursion rela¬ 
tions applied at the level of Feynman diagrams. In particular, one loop diagrams are written in 
terms of tensor integrals multiplied by coefficients that are determined numerically. The tensor 


integrals are then reduced to scalar 
implements the algorithms of Refs. 


i nteg rals by means of the COLLIER library 


151 


150] which 


152]. Rational terms and renormalization countert¬ 


erms are included as additional tree-level Feynman rules. MUNICH and SHERPA compute 
real QCD/QED correction s and pe rform the subtraction of infrared singularities within the 


dipole sub 
In Ref. 


raction scheme 


153.1154]. 


137] the one loop corrections of absolute order a 2 s a 3 to Z +2 jets have been computed 


using the code RECOLA. Also in RECOLA the computation of LO and NLO amplitudes is 
based on recursion relations which however act on off-shell currents (that correspond to sum of 
Feynman sub- diagrams), leading to a different treatment of color with respect to OpenLoops. 
As in OpenLoops, the algorithm constructs numerically in four dimensions the coefficients of 
the tensor integrals that are then reduced using the COLLIER library. Rational terms (as well 
as UV renormalization counterterms) are introduced at the level of LO amplitude generation, 
while IR singularities are canceled in the dipole scheme. 

One loop electroweak and QCD corrections to the processes pp —>■ it + H , pp —>■ tt + 
W an d pp —> i t + Z have been computed in the MadGraph5_aMC@NLO framework in 


Refs. 


m 


in MADFKS 


139 [. IR singularities are subtracted in the FKS 


155 


156] scheme implemented 


157 ]. while one loop amplitudes are generated by MADLOOP which provides 


an independent implementation of the OpenLoops algorithm (where one loop integrals are 
computed either via tensor integral reduction or by means of the OPP method). 

As a final remark, it is worth noting that the automation of one loop EW corrections is still 
in its early stages, but rapidly evolving. 


II. VIRTUAL CORRECTIONS AND SCHEME DEPENDENCE 


In one formula, the basic ingredients of an NLO computation can be sketched as 

_ _Virt. CT , T 
(7„ — <J n + (J„ + (X 


NLO _ ^Virt. | ^CT | ^.Real, sub > I PRrad / Real sub 

/ blAU lC7 1 1 Vn 


'n+1 u n+1) "> 


(i) 


where n is the number of external legs of the LO process under consideration. In Eq. © , ovf irt ' 
is the contribution of the diagrams with one virtual loop: these diagrams can be ultraviolet 
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(UV) divergent when the number of propagators involved in the loop integral is smaller than 
four, so that a regularization procedure is needed in order to keep track of the UV singularities 
that will be removed by means of renormalization (in additive renormalization prescriptions, 
this is achieved by adding the appropriate counterterms -cr^ T in Eq. ([I])- to the virtual correc¬ 
tions). Besides the UV singularities, in massless gauge theories such as QED or QCD, cqf irt ' 
can be also infrared (IR) divergent: the 1R singularities of the virtual one loop corrections are 
cancelled by the corresponding IR singularities that appear in the computation of real NLO 
corrections. In Eq. ([I]) we assume that the integral over the real phase space has b een regular¬ 


ized by means of a subtrac tion procedure |158l |. such as the Catani-Seymour |l53l . 154 ] or the 

156], and crR eal > sub represents the integrated contribution of the singular 


FKS subtraction 


155 


part of the real NLO corrections. In the context of NLO corrections to hadronic processes, 
additional singularities arise due to the collinear splitting of massless initial state partons: as 
these initial state singularities are universal, they can be absorbed in the parton distribution 
functions (PDFs) and have to be consistently subtracted form the one loop calculation by 
adding to Eq. (JT]) the corresponding collinear counterterms. 

The commonly used regularization prescription for UV singularities is dimensional regular¬ 
ization (DR) 159] . The basic idea of DR is to promote the number of space-time dimensions 
from 4 to D = 4 — 2euv < 4, so that the loop diagrams read 

(2 ^y~ D 


d N(q, {Pi,^(Pi)}) 


m* 


d q 


n;io [(g-u) 2 -m?r 


( 2 ) 


where q is the loop momentum, {pi,ip(pi)} schematically denote the external momenta (con¬ 
sidered as incoming) and their spinors or polarization vectors, rq = JV =nPii m i * s ^ ie mass 
of the particle in Z-th loop propagator and /i is a mass parameter introduced in order to keep 
the dimension of the integral fixed for varying D. In DR singularities appear as 1/euv poles 
in the limit D —> 4. In Eq. ([2]) the full Lorentz and Dirac structure of the loop diagram is 
included in the numerator AT. In QCD calculations, DR is the standard regularization pre¬ 
scription also for IR singularities (even if D = 4 — 2ei R > 4 and e IR = —euv). In the context 
of NLO clectroweak computations, the traditional IR regularization scheme is mass regular¬ 
ization (i.e. IR singularities are regularized by giving an unphysical mass to the photon and 
the massless fermions), however DR is becoming the new standard procedure in order to treat 
the IR limit of QED and QCD on the same footing (this is of particular importance not only 
for the computation of mixed QCD and electroweak corrections, but also for the extension to 
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clectroweak computations of NLO +pa rton shower matching algorithms originally developed for 


QCD calculations |33|, 


36 


59 


61 


160j]) 


In the original formulation of DR in Ref. [1591 ]. not only the loop momentum but also all the 
other momenta become D-dimensional together with the Lorentz covariants (such as and 
y M ), which become formal objects obeying specific algebraic relations. This kind of continuation 
to D dimension is consistent and well defined for non-chiral theories, while problems arise in 
the t reat ment of the y 5 matrix and the £^ pa antisymmetric tensor. More precisely, already in 


Ref. 


159 ] has been pointed out that the main properties of the y 5 matrix in four dimensions (i.e. 


{y 5 , 7 m } = 0, Tr(y 5 y /i y i/ y p y <T ) = 4 i£ pvpa together with the trace cyclicity), in D dimension could 
lead to algebraic inconsistencies in presence of anomalies. In the literature, several different 
recipes to handle the y 5 matrix in the context of DR have been proposed: 

• ’t Ho oft y 5 prescription. This y 5 prescription has been proposed by ’t Hooft and Velt- 


man 


Maison 


1591, A kyeampong and Delbourgo jl6ll . 162] and systematized by Breitenlohner and 


163]. According to Ref. [1591 ]: 


7 5 = hVtV = 


(3) 

{7 5 5 7 M }^0, {y 5 ,y^}=0, [y 5 ,y^]=0. (4) 

In equation (j4j) we introduced the notation y M = y^ + y^, where y^ and y M represent 
4 and D — 4 dimensional objects, respectively. Even if algebraically consistent, the ’t 
Hooft y 5 prescription leads to spurious anomalies that should be removed by introducing 
appropriate finite counterterms. 


KKS prescription. The scheme proposed in Refs. 


164 


166] preserves the anticommutation 


relationship {y 5 , y^ 1 } = 0 and prevented cyclic property in Dirac traces to avoid algebraic 
inconsistency. A projection on four-dimensional subspaces is needed in the redefinition 
of trace operation: 


TtnVi'W) = r tt(j , (7 5 b'‘7‘W), p<7 5 ) = jfw rrrr- (5) 


Since no cyclic permutation of Dirac matrix inside traces is allowed, all the traces must 
be read starting from the same vertex (reading point prescription). 
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Naive DR. In this prescription { 7 s , 7 M } = 0, Tr('y 5 'y tl 'y l, 'y p 'y (T ) = and cyclic per¬ 

mutation of gamma matrices inside traces are allowed (see for instance Ref. (l67| ). Even 


though these properties can be incompatible, the resulting ambiguities are expected to 


cancel in non-anomalous theories 
coimterterms from Ref. 


168 


170]. We will take the formulae for the electroweak 


171], where the naive prescription has been used: accordingly we 


will use this prescription in the examples of subsection III Dl for our CDR (conventional 
DR) results. 

These three different prescriptions for the treatment of 7 s are three possible examples of different 
realizations of DR: all of them share the same definition of the denominator of the loop diagram 
in Eq. ([2]), while they differ in the treatment of the numerator of Eq. ([2]). Different choices about 
the numerator algebra in Eq. ((2]) concerning the dimension of the momenta, the Dirac matrices 


or the spinors and polarization vectors, c 
conventional DR scheme (CDR) 


159 


172 


ehne different DR schemes, such as for example the 


(HV) scheme [159] or t 


174] (or some variant thereof), th e’t Hooft Veltman 


re Four Dimensional Helicity scheme (FDH) |175I . 176]. Dimensional 


Reduction (DRED) |177l -180| can be regarded as a DR scheme in which all the objects appearing 


in the numerator of Eq. ([2]) are four dimensional (except, of course, the loop momentum): this 
is the regularization scheme used by GoSam for the computation of one loop diagrams and 
it is particularly appealing because it allows to work out the numerator of Eq. (J2]) using four 
dimensional relations plus some trivial algebra to treat the loop momentum components in the 
extra dimensions without introducing special Feynman rules (for the hnite integrals involved 
see for example Ref. M). 

Different DR schemes lead to different results for the computation of one loop diagrams. The 
differences consist of e-dependent terms, that vanish in the D —> 4 limit unless they multiply the 
1/e poles, giving rise to hnite terms (usually called rational terms in the literature). The main 
purpose of the present paper is to show how the automation of NLO electroweak corrections 
can be achieved starting from existing tools, such as GoSam (for the computation of one 
loop matrix elements) and MadDipole (for the IR subtraction). If we recall Eq. (JT]), one 
of the technical issues we have to face is the matching between the cqf irt ' term computed by 
GoSam in DRED and the contributions coming from renormalization and from the integrated 
dipoles (<t^ t and cr^ eal,sub , respectively) that are usually computed in the CDR scheme in the 
literature. Although the different building blocks of an NLO calculation are regularization 
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scheme dependent, as pointed out in Ref. 1811 ] unitarity implies that the regularization scheme 


dependence vanishes in the sum of , a^ T and cr^- eal ’ sub , once these three terms have been 
computed within the same regular izati on scheme. For the counterterms part, in particular, we 


implemented the formulae of Ref. 


17l| for the on-shell renormalization scheme and derived the 


conversion factors from CDR to DRED. The infrared subtraction procedure implemented in 


MadDipole is based on the Catani-Seymour dipole formalism: 
factors for the integrated dipoles have been taken from Ref. 


he CDR to DRED conversion 


181 ] (with a straightforward 


generalization to the case of QED radiation). In subsection III Al we give a short overview 
about the electroweak renormalization at one loop in the on-shell scheme, we also explain our 
treatment of massless fermions and list the relevant regularization scheme conversion factors for 
the electroweak counterterms. The subtraction of IR singularities is described in subsection III Bl 
while in subsection III Cl we address the issue of the regularization scheme dependence of the 
infrared subtraction procedure. As examples on how the things work, in subsection III Dl we 
show the CDR and DRED results for the NLO electroweak corrections to three simple processes: 
Z —> uu, W + —> ud and uu —> cc. In particular, with the examples we show the regularization 
scheme independence of the sum of virtual one loop corrections, renormalization counterterms 
and integrated dipoles when these three contributions are consistently computed in the CDR 
or DRED scheme. 


Electroweak renormalization and counterterms 


We use the on-shell scheme 


182l - ll9l| for the electro wea 


the sake of simplicity, we employ the notation of Ref. 171j and refer to Appendix B of 171 ] 


c renormalization at one loop. For 


for the explicit expressions of the self energies involved in the formulas for the electroweak 
counterterms. Each free parameter of the Standard Model Lagrangian is considered as a bare 
parameter and it is split into its renormalized version and the corresponding counterterm, 
namely: 


e 0 



e + be = (1 + 5Z e )e, rrif : o = rrif + 8m f 

Ml, + 6 Ml. M| 0 = M| + 6 Ml M 2 H0 = M 2 „ + 6M 2 H , 


( 6 ) 
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(where / stands for a generic fermion field), then for each external field we introduce the 
corresponding wave function renormalization counterterm: 

Wo = (1 + \8Z w )W, H 0 = (1 + Uz h )H, 

f 0 L = (1 + \sz))f L , f 0 R = (1 + lszf)f R , 

Z 0 — (1 + -5Zzz)Z + —SZzaA, 

Ao = (1 + — 5Zaa)A + -8ZazZ. (7) 

This way the Lagrangian can be split into the basic Lagrangian and the counterterm Lagrangian 
Co = C + SC, where SC provides the Feynman rules for the counterterms. 

In the on-shell scheme the form of the renormalization counterterms is fixed by the following 
requirements: 

• the renormalized masses are the real parts of the poles of the propagators at one loop; 

• the residues of the renormalized propagators are equal to one; 

• the renormalized electric charge is obtained from the eey vertex in the Thomson limit; 

• the renormalized Higgs tadpole is set to zero. 

With the above conditions the one loop EW counterterms for fermions, gauge bosons and EW 
couplings read: 


Sm, = ^Re (s{; L (mj) + E/;" (m'j) + 2 T,{f Imjj 'j , 

/>z!; L = -ReEV (mj) - mj-C (r eS {: L ( p 2 ) + ReE^ R (p 2 ) + 2ReEf (p 2 )) 

SZ{: R = -ReE^ s (mj) - mj-C (ReEV (p 2 ) + ReE^ (p 2 ) + 2ReEf (p 2 )) 


( 8 ) 
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SM 2 W = ReS^T (M 2 ,) , 

SM 2 Z = ReS| z (M 2 Z ) , 
d 


SZ\y — 


SZzz — 


d 


ReE^ (/c 2 ) 

S MM Z (* 2 ) 


dk 2 

d 


k 2 =M 2 , 


k 2 =Mg 


5Z aa = 


dk 2 


SZ AZ = —2Re 


S AA (fc 2 ) 

S? z (Ml 


fc 2 =0 


SM 2 H = ReZ H (Af£) , 


= 2 


= 


E AZ (0) 

Ml ’ 


5 


<% 2 


ReE H (/c 2 ) 


k 2 =Mjj 


(9) 


^ C 7 % 1 r 7 1 & x^AA (l 2\ 

— = ~A SZ AA - -SZ ZA = -— -E r (k ) 

e 2 cw 2 ° Ah2 v 7 


2dk 2 ^ T 

Sc w i ( m 2 w _ sm 2 z \ = i ~ (y% (m 2 ,) 
cw 2 V M 2 , Mf y 2 e 

^SlV C^r 


k 2 =0 

\pZZ 

Z-Jjn 


s w S AZ (0) 
Cw Mf 


(m: 




M 2 


>w 

sw sfy c vk 
5t = -T h . 


1 c 2 ,„~ (T^(M 2 


rw 
9 s 2 


Re 


wv 


MJ, 


s| z (mi; 

Mf 


( 10 ) 


Where the notation Re means that the real part is taken only for the scalar functions con¬ 
tained in the unrenormalized self energies. As in our computation the CKM matrix Vqkm is 
chosen to be diagonal and no renormalization of Vckm is needed, we do not discuss here the 
renormalization conditions for the quark mixing matrix. 

In the following we list the conversion rules from the CDR to the DRED regularization 
schemes. For electroweak couplings and gauge bosons the conversion rules read: 


SZ) RED 

$ z aa ED 

SM^ red 

SZ% RED 


CDR 




&Zaa R 


a 1 
47T 3’ 
a 2 
Rr 3 ’ 


SZ 


ZA 


= SZ 


ZA 


SZf^ ED = SZf) R 


oi 2 c 2 v 
4vr3s^’ 


a 2 M^r 


r M 2 CDR 

01V1 W 2 

47t3 s z 


w 


SZ 


DRED c yCDR , a ^ C W 


AZ 


= SZ 


AZ 


47t 3 sw 


SZ ? DR 


SMl DRED = SM 2 7 cdr + 


Ml 


J H 


SM 2DRED = 5M 2CDR + “ 3 ( 2c w™w 


SZ ERED = SZ IP R 

a 2 M 'Dy 
4vr3 s 2 w 
+ Ml) 


a 


4tt3 4 / 


47T 


2r 2 s 2 


( 11 ) 
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The CDR to DRED conversion rules for up-type quarks, down-type quarks and leptons, re¬ 
spectively, are: 


k DR£D 

DRED 

j^RDRED 


<5m u CDR + 
5Z RCR>R - 

SZ RCDR_ 


a (3 — 4 s^y) 2 + 2c^(9 + 8sj+) 

72^ 

ot (3 — 4 s^y) 2 + 2c\y(9 + 8s^/) 
4tt 36c^s^ 

cc 4 
4tt 9c^ ’ 


( 12 ) 


5m d DRED 


SZ L DRED 


5Z r dRed 


5m d CDR + 


5Z rcdr 


SZ RCBR 


a 9 + 4s^ + 2c^(9 + 25^+) 
4vr 72c w s w 

a (3 — 2 s^y) 2 + 2cj+(9 + 2sj+) 
4vr 36cj+s^ 

a 1 
4tt 9cj+ ’ 


(13) 


5m^ u = 5m e 


CDR 


a 1 + 4 s^/(s^/ — 2) + cfyl 2 + 4 s^) 


+ m - 
47r 


8 o 2 s 2 


DRED _ SZ^ CDR 


^^RDRED _ SZ^ CDR 

figL DRED _ ^LCDR 


a 1 + 2c 2 


w 


4vr 4c^s^ 
a 1 
4vr c 2 ,’ 


a 1 + 2c 2 


w 


4tt 4c^s^ 


(14) 


As in our numerical examples we deal with external massless quarks, following for example 
Ref. [59}, we adopt a slightly modihed version of the on-shell scheme where the light quark 


masses are set to zero in the corresponding external wave function renormalization counterterms 
(and mass counterterms for light quarks vanish). In this approach eqs. ffl2|) - f|T3]l for massless 
up and down-type quarks in particular become: 


SZ LDRED 

= SZ RCBR - 

a (3 — 

4s 2 ,) 2 + 18c4 

47T 

36c^s^ 

fo Z R DRED 

= SZ RCDR - 

oi 4sj+ 
A-k9c 2 w 

J 

SZ L DRED 

= 8Z L d CRR - 

a (3 — 

2S^/) 2 + 18 Cyy 

47T 

36c^s^ 

6Z RVRED 

= 8Z rcdr - 

a s 2 w 


4vr9c^ 



(15) 


(16) 
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Large logarithms of the fermion masses appear in the counterterms SZ e and SZaa■ They are 
basically related to the running of a from the Thomson limit to the electroweak scale. These 
terms are singular in the limit m q —> 0, however they can be systematically cancelled by 
following the recommendation of the EW dictionary of Ref. [62], namely by choosing a(0) as 
coupling constant for the LO vertices involving external photons and a(M|) or oiq for the 
remaining LO vertices. In the former case the logarithms of the light fermion masses cancel 
in the combination SZ e — ^SZaa, while in the latter they vanish as a re sult of the finite shifts 


SZ e \ a (M 2 z ) £Z e |a(o) — TjAa(Mf), 8Z e \ G/i —> 5Z e | a ( 0 ) 


iAr 


49, 


119 


17 If . The remaining 


dependence on the light quark masses in the counterterms vanishes in the m q —> 0 limit. 
Indeed, in this approach light quark masses can be safely neglected from the very beginning. 

We conclude this section with a short remark on unstable particles. In general the description 
of resonances in perturbation theory requires a Dyson summation of the self energy insertions. 
This leads to a mixing of perturbative orders that could break gauge invariance. In the context 
on one loop electroweak corrections unstable particles are usually treated within the so-called 
complex mass scheme (CMS) 


192 


194 ], In the CMS unstable particles masses are promoted 


to complex numbers through the replacement 


rriy —> /ay = niy — imyY v , 


(17) 


(where Ty is the decay width of the unstable particle V) and the complex Weinberg angle is 
then defined via 

cos 2 d w = . (18) 

hz 

Within the CMS renormalization can be performed in a modified version of the on shell scheme. 
The renormalized masses of the unstable particles are defined as the poles of the corresponding 
propagators in the complex plane. The resulting expressions for the counterterms are basically 
the same as in eqs. (IHl fTOj) with the important difference that all the self energies involved 
become complex (via their dependence on the complex couplings and masses) and no real part 
is taken. In the CMS also the momentum flowing in the self energies is a complex variable. 
However, it is possible to realize a minimal version of the on shell renormalization in the CMS 
by expanding the self energies about the real arguments in such a way that one loop accuracy 
is retained. This in particular requires the introduction of an additional renormalization factor 
for the W (and also the top) self energy in order to deal with the branch cut at k 2 = ^ 
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arising from the photonic loop diagrams. The resulting counterterms are still complex, but 
depend on real momenta. 


B. Electroweak insertion operator 


In equation ([T]) we briefly sketched the basic building blocks of an NLO calculation. We 
pointed out that virtual corrections cqf irt ' and real contributions are separately infrared 
divergent and the IR singularities only cancel in the sum of the two terms (as far as infrared 
safe observables are considered and the initial state collinear singularities are properly absorbed 
in the PDFs in the case of hadronic processes). In Eq. ([I]) we also introduced the infrared 
subtraction term cr^’, and its integrated counterpart orR ea l sub . a^ 1 is a function of the n + 1 
kinematics that behaves as the matrix element for the real corrections in the soft/collinear limit: 
this way — c r ^ b i in Eq. (jT)) is IR finite and can be integrated numerically, should 

also have a simple expression, so that the IR subtraction term can be integrated analytically 
over the real radiation phase space: the resulting integral er^ ca1, sub is a function of the tree level 
n particle kinematics that has to be added back to the virtual corrections in order to cancel 
the IR poles of the one loop contributions. 

Owing to the universality of the infrared limit of massless gauge theories like QCD it is 
possible to develop process independent algorithms to build the IR subtraction terms and 
of course their integrated counterparts (the dipole subtraction or the FKS methods are two 
examples of such algorithms). The infrared properties of one loop amplitudes are well 


Known 


and in the context of QCD are usually described by a so called insertion operator 


195 ], that 


corresponds to the integrated subtraction terms described above. These considerations can 
be easily translated to electroweak one loop amplitudes and the infrared properties of such 
amplitudes are then described by an electroweak insertion operator. The electroweak infrared 
insertion operator is particularly useful as it allows to check the cancellation of double and 
single pole of the renormalized ampli tude for a given phase space point. It can be obtained 


for instance from the results of 


153 . 


154 ] by, simply speaking, replacing the color factors by 


electric charges. More specifically we start by dehning the singular pieces as: 

k ' sing ,Q E ik i ^Ifc) T ~k^ e (Sjfc, UTj, TTlk} \ ' l-AT^ornl • (19) 


The index i denotes the combination of emitter and unresolved particle, the index k denotes 
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the spectator. The coefficients of the double and the single pole depend on the masses and are 
given by 


rrii > 0 , m k > 0 : 


V t 2 (s ik ,mi,m k ) = 0 


V € (s ik ,mi,m k ) 


logp y/X (si k ,m%,m%) 1-Vik 

-> V ik = -2-2 ’ P = \ Tb-r» 

Vik Sik -mf-ml V 1 + Vik 

A(x, y, z) — x 2 + y 2 + z 2 — 2 xy — 2 yz — 2 xz, 


mi > 0 , m k = 0 : 

V e 2(s ik ,mi,0) = i 

1 f 777 ? \ 

V t (s ikl m t , °) = - y- + log —±J 

mj. = 0 , m k > 0 : 

V e 2(s ik ,0,m k ) = ^ 

V £ (s ik , 0, m k ) = \(^ + lo S 

mi = 0 , m k = 0 : 

V e 2 (s ifc , 0 , 0 ) = 1 

V e (s ifc ,0,0) = ^. (20) 

The insertion operator depends on the fact whether emitter/spectator are in the initial or final 
state. In particular, we have: 


initial-initial: 


j / \ _ ^c^iQi^kQk / \ | i i |2 

le 2 \$iki TTlii n h e 2 \Siki ITlij 1Tl k j ' - / 'dj3 orrl 

T ( 2\ n c&iQi&kQk^, / \ (3 . h \ | , j 12 

I e {s ik ,mi,m k ,n ) =- V e {s ik ,mi,m k ) - + log— • \M B om\ 

e \2 s ik ' 


initial-final / final-initial / final-final : 

T / \ _ ^c^iQi^kQk / \ | . , | 

J-c 2 \&iki ^77-i? ^72fc) 9 * 


I e (s ife ,m i ,ra fe ,/h 2 ) = n<(Tl(Il(Tk(lk (v e (s ik , mt, m k ) + log—V e 2 (s ik , mi, m k ) 

^ V 


• |Af 


Born 


(21) 
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The factors g*, g*, denote the fractional electric charges of emitter and spectator and n c is the 


num 


jer of colors that can occur in the splitting (i.e. either one or three). In addition, we follow 


Ref. 


196] and introduce the sign factors cy*. which are defined to be +1 for incoming fermions 


and outgoing anti-fermions and —1 for incoming anti-fermions and outgoing fermions. This can 
be extended to the case of a final state bb-boson and yields a +1 for an external W~ and a 
— 1 for an external W + . The sign factors ensure the conservation of the electric charge which 
is given by 

<?n ■ Qn = 0 . ( 22 ) 


C. Subtraction terms and scheme dependence 


Numerical simulations require the use of subtraction terms to render the real emission con¬ 
tribution finite. One of the widely use d subtraction methods in the context of electroweak 


calcu 


1961 . 


ation is the dipole formalism 


153] which has been adapted to electroweak calculations 


197 ], In practice the subtraction methods usually also takes care of issues regarding the 


regularization scheme by adding appropriate additional terms to cancel the scheme dependence 
in the virtual corrections. 

As mentioned above, different regularization schemes lead to differences in the finite con¬ 
tribution for the virtual corrections. At the one-loop level the transition from one scheme to 
the other can be obtained by a simple shift that is proportional to the Born contribution 181]. 


Following the derivation of Ref. 
term of the form: 


1811 ] . in the case of photon emission one obtains a conversion 
1 


$RS = 


iQi&i Qk Ok 


(23) 


changing from CDR to DRED. Here i denotes the emitter, k denotes the spectator, g^*, are the 
fractional charges of emitter and spectator respectively, and oy^ are the sign factors as defined 
in Eq. (1221) . However there is only a contribution in case of massless emitters. Soft singularities 
are independent of the scheme and therefore a scheme dependence only arises from collinear 
splittings. As these are finite in the massive case there is no transition term for massive emitters. 
Eventually the full transition term is obtained by summing over all possible emitter-spectator 
pairs. 
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FIG. 1: One loop electroweak corrections to V —» f f. The diagram on the left gives rise to the V a 
contributions, while the V& terms come from the diagram on the right. 

D. Examples 

1. Z —» vv 


As a first simple example of the applications of the formulas of the previous subsections we 
consider the decay Z —y vv. At the LO the corresponding matrix element is: 


M 0 = ieg u u(pi)'y fl uj-v(p2)£f,(pi +p 2 ), 


(24) 


where gj = (// — s^Q /)/ swcw, w± = (1 ± y 5 )/2 and is the polarization vector for the 
on-shell external Z boson. 

The two classes of diagrams contributing to the one loop electroweak corrections to Z —y vv 
are shown in Fig. [Q The unrenormalized O(a) corrections to the Z decay in neutrinos read: 


Q, f 

A4 V irt. = Mo— s 4 a 1 2 V a (0,M|,0,M z ,0,0) 

4tt l 

+ sJ ^ V a (0, M|, 0, M w , 0, 0) 

+ J V b (0,Ml0,0,M w ,M w )}, (25) 


where we used the notation of Ref. 


171] for the vertex functions V a {b)- It is worth noting 


that Eq. (l25jl does not depend on the regularization scheme. What is scheme dependent is the 
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expression of the V a (b) functions. When performing the calculation in CDR we find: 


V. 


GDR (0,s,0,m,0,0) = -2sC 0 (0, s, 0, m, 0, 0) (l + — J - B 0 (s, 0, 0) (3 + 


.77T 


+2B 0 (0,m,0)(2 + 


m 


- 2 


(26) 


2 2 
m\m 2 


V 6 CDR (0,s,0,0,mi,m 2 ) = 2(Vrq + mj + ) C 0 (0, s, 0, 0, m 2 ) 


-B 0 (s,m 1 ,m 2 )(l + 


m\ + m\ 


+ ^2 H—-iji?o(0, 0 , 7777 ) + ^2 + —^B 0 (0, 0,ml), (27) 

(where Bo and Cq are the usual Passarino-Veltman two and three point functions, respec¬ 


tively 


198 


200j), while the DRED results are related to CDR ones as follows: 
V a DRED (0, s, 0, m, 0, 0) = V a CDR (0, s, 0, m, 0, 0) + 1 
V fe DRED (0, s, 0, 0, mi, m 2 ) = V GDR (0, s, 0, 0, m l7 m 2 ) + 1. 
The O(a) corrections to Z —y vv coining from renormalization read: 


M C T = Mo{bz z + 5Z!; + 5Z e + ( s -f- 


(5s 


w 


Sw 


(28) 


(29) 


Using Eq. (125]) and the scheme conversion rules for the counterterms in eqs. dill) and (ITT]) . it 
is easy to show that ( - Mvi ^ + Q - Mc!T ) CDR = ( ^ V1 % + 0 ^ CT ) DRED - 

Note that Eq. (1291) contains large logarithms of of the light fermion masses, that can be 
removed by using ck(M§) as input parameter in Eq. (1241) and performing the corresponding 
finite renormalization SZ e \ a {M%) ^^ e U(o) — |Ack(M|), where 


y\AA /Q 2 \ 

A Q (Q 2 ) = n^ %h ,(0)-Re(nyy i , lt (Q 2 )), n^, i8hl (Q 2 ) = /£ “g‘ , (30) 

and ^/e light represents the contribution of light fermion loops to the photon self energy (which 
has the same expression both in DRED and in CDR). Once the singular terms in the m q —> 0 
limit have been subtracted, the remaining dependence on the light quark masses can be safely 
neglected. 


2. W + ->• ud 


With the notation of Eq. (|24p . the matrix element for the process W + —> ud is: 


Mo = i 


y/2s\y 


u(p 1 )'fM-V{P2)s f i(Pl +P2)- 


(31) 
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The one loop electroweak corrections to W + —> ud, are given by the two classes of diagrams 
in Fig. m At variance with the Z —> vv case, in this case also photon loops contribute to the 
corrections. We basically follow the computation of Ref. 1711 ]. with the important difference 
that we use dimensional regularization for the IR singularities coming from photon loops. These 
1R singularities are cancelled by the QED dipole subtraction terms described in subsection IIIB1 
The unrenormalized virtual one loop matrix element is: 


Afyirt. = A4o—| QuQd V a (0,M^,0,0,0,0) 

+9 U 9d Mwi 0 ) M z , 0 , 0 ) 

+Qu V b (0,M^,0,0,0,M w ) 

—Qd V b (0,M^,0,0,M w ,0) 

+9uV b (0,M^0AM z ,M w ) 

~9dZ V b {0,M^,0,0,M w ,M z )Y (32) 


while the contribution coming from renormalization can be written as: 

Mct = M 0 \\5Z W + Uz^ + Uz T d +SZ e -^\. 

1/ Z Z S\y J 


(33) 


As in the previous example, eqs. (132113 3 p are completely general, while the expression of the 
vertex functions and the counterterms depend on the regularization scheme. If we define h Virt ' 
via the relation: 


2 Re 


A4vi rt . + -Mct ) Afo 


_ a jVirt. 

27T 


I M c 


(34) 


from eqs. (J28J), (I15H16() and (HIT) , we find that the one loop virtual electroweak corrections have 
the following regularization scheme dependence: 


rVirt. rVirt. _ 

°DRED °CDR — 


\qi+\qi ■ 


(35) 


Also in this case, the scheme dependence of the one loop electroweak corrections is canceled 
by the corresponding scheme dependence of the QED integrated dipoles that are obtained 
from Eq. (f23l) by summing over all emitter and spectator pairs. As a result, the sum of the 
renormalized one loop virtual corrections and the infrared subtraction terms is regularization 
scheme independent. As in the previous example, the dependence on the light quark masses in 
Eq. (l33ji can be removed by using a(M z ) or cxg m as input parameters. 
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FIG. 2: LO and virtual NLO electroweak diagrams for the process uu —>• cc. 

3. uu —>• cc 

As a last simple example we consider the one loop electroweak corrections to the O(as) 
process uu —* cc. As the 0(a) corrections are hclicity-dependent, it is convenient to write the 
LO matrix element as: 

Mo = ( 9s * 3 k * u(Pc)7 l *wxfV(p s )v(pii)'y l j.wx i u(p u )= M>(Ai,A/). (36) 

KPu+Pu) KXf=± Aj,A/=± 

Fig.[2]represents the diagrams contributing to the NLO electroweak corrections to the process 
uu —y cc. As in the previous example, we regularize dimensionally the QED-like IR singularities 
coming from the photon loops and also the QCD-like IR singularities associated with the gluon 
exchange between the fermionic currents in the initial and in the final states. For a given 
helicity configuration, the virtual electroweak corrections read: 

M Virt .(\,A f ) = M 0 (\,A f )^{s Vertex (\) + S Vertex (Af) + <5 Box (A,, A/) j, (37) 

with 

£ Vertex (a) j (^ )2 ^ + ^ Qj ^ Q, Q) 

+QI v a (0,(p u + p s ) 2 ,0,0,0,0) 

V a (0,(p u +p s ) 2 ,0,M w ,0,0)}. (38) 

The counterterm contribution is: 


(39) 
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where 5Z+ = 5Z * <5Z" = and 6Zu {R) = 5Zc {R) . 

It is worth noting that the box contribution in Eq. (157]) does not depend on the regularization 
scheme (CDR/DRED) and d oes n ot give rise to rational terms neither of UV type (since it 
is UV finite) nor of IR type [201 ]. As a consequence, using eqs. (125]) in (155]) together with 
the scheme conversion formulas for the counterterms C0 we can show that for each of the 
helicity configurations the one loop electroweak corrections depend on the regularization scheme 
according to: 

^DREd(^*> ~ ^CDR.(^*> ^/) = ^Qui (40) 

where d Virt '(A ? ;, A/) corresponds to Eq. (155]) with the replacement A4 —> M.(Xi,Xf ) (of course 
Eq. <m a ^ so holds at the level of unpolarized matrix elements). Also in this case, the scheme 
dependence of the one loop electroweak corrections is cancelled by the corresponding scheme 
dependence of the QED integrated dipoles. 


III. REAL EMISSION AND SUBTRACTION PROCEDURES 


Conceptually the calculation of the real emission processes and the appropriate subtraction 
terms for the cancellation of QED infrared singularities in numerical simulations is very similar 
to the case of QCD calculations. Many concepts and methods that were originally developed for 
QCD calculations can easily be adapted to electroweak computations. The additional radiation 
of a photon is a source of infrared singularities in regions of phase space where the photon 
becomes either soft or collinear to other charged massless particles. For sufficiently inclusive 


observables, the KLN theorem 


202 


203] guarantees the cancellation of infrared singularities 


from final state particles against their counterparts from the virtual contributions. Collinear 
initial state singularities have to be either absorbed into the parton distribution functions, or 
these singularities have to be regulated using the physical mass of the incoming particles. The 
corresponding logarithms of the masses are a potential source of large corrections. 

In the context of one loop electroweak calculations, the real contributions are basically 
defined as the real QED corrections to the process under consideration. Nevertheless, also 
the radiation of an additional W/Z boson would contribute at the same order in perturbation 
theory. In the literature, the extra emission of additional W and Z bosons usually is not 
included in the computation of real O(a) corrections for two reasons. First of all the gauge 
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bosons decay, so that in principle they lead to final states which are not degenerate with the LO 
ones. The second reason is that real purely weak corrections are always finite: in fact, even in 
the Sudakov regime, where the corrections are dominated by the Sudakov logarithms which are 
the IR limit of virtual weak corrections, the gauge boson masses act as physical IR regulators, 
ffowever, in those kinematical regions where the Sudakov corrections become large, the effect 


of the partially compensa 


contributions. In Refs. 


90 


ing radiation _of real gaug e bosons may lead to significant positive 


92 


24 


135 


139 


204 


208] the issue of real weak boson emission has 


been addressed in analogy with QED or QCD: all the diagrams obtained from the LO ones with 
the emission of an additional W or Z boson are considered as real corrections, the additional 
gauge boson is produced on-shell and integrated over the full phase space. The result is a 
significant cancellation between real and virtual corrections, which however may be incomplete 


due to the incomple 
violating effects [90 


e average on the isospin of the initial state particles 


243] (Bloch Nordsieck 


92[). Real weak boson emission processes have been considered from a 
more phenomenological point of view in Refs. R Lot! l209i 21ol ] : the additional gauge bosons 
decay and they are included in the real corrections only when the final states are degenerate 
with the LO ones. Also in this case, the cancellation between real and virtual weak corrections 
is only partial and strongly dependent on the event selection under consideration. Recently, two 
different implementations of multiple wea k boson emission processes in the context of parton 


showers have been presented in Refs. 


2111 ] and [2121 ] . 


From the theoretical point of view the real emission of a photon can be handled in the 
same way as it is done in QCD, namely by defining jets as infrared safe observables, and a jet 
may or may not contain a photon. A clean way of doing that is the ’democratic clustering’ 
approach [213j where all the partonic final states (hadrons and photons) are treated on equal 
footing. In that sense the photon is then treated fully inclusively and the results are infrared 
safe. For many final states however the presence of an additional photon can experimentally 
be tagged and therefore be distinguished from the original Born process. For instance a muon 
being collincar to a photon can be disentangled as they end up in different parts of the detector. 
From the theoretical point of view, this implies that the photon is not treated fully inclusively 
but parts of the collinear singularity is cut away. This leads to the problem that the collinear 
singularity in the virtual corrections is not fully canceled by the corresponding one in the real 
emission. As the muon mass acts as physical cut-off for the collinear singularities, the only effect 
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of the miscancellation is a large QED correction. However, the incomplete cancellation would 
be an issue in the case of the photon radiation off a massless particle. There are two options to 
handle this problem. The first starts from the observation that a photon can either be produced 
in the hard scattering, i.e. it can be treated perturbatively, or it occurs much later at lower 
energies, when hadrons are created and decay into more stable hadrons. This happens in an 
energy re gime that is non-perturbative and which can be described by the photon fragmentation 


function 


213], The photon fragmentation function is an observable that has to be measured 


experimentally, similarly to the parton distribution functions. And in the same way as initial 
state collinear singularities can be absorbed into the parton distribution functions, the part of 
the collinear singularity that is not canceled by the virtual corrections can be absorbed into the 
photon fragmentation function. The second option to deal with this incomplete cancellation is 
to (also experimentally) only cut away finite pieces such that the cancellation of the singularity 
remains unspoiled. This can be achieved by reducing the allowed amount of photonic energy 
within the jet as a function of the R-separation between photon and tagged particle such that 
the photon energy vanishes in the limit of collinearity 2141 ]. We note that this problems also 
occur in fixed order QCD calculations where photons are present. 

Additional complications arise from the fact that for a general process one cannot distinguish 
between QCD corrections and electroweak corrections but all contributions at a given order in 
a and a s have to be taken into account. For the real emission this includes contributions from 
a gluon exchange between the interference of a QCD Born diagram with an electroweak Born 
diagram. It is worth to stress that these contributions are therefore not the absolute squared of 
an amplitude but just the interference term between two different types of amplitudes. They 
have different properties concerning their color -but also their kinematical- structure compared 
to standard tree-level like contributions. These terms can be handled with standard methods 
from QCD and we refer to section [TV] for a concrete example, nevertheless their automated 
generation is in general not trivial (but it is within the reach of the tools dealing with automated 
calculations). 
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IV. ELECTROWEAK CORRECTIONS TO W + PLUS TWO JETS 


In the last part of this review we show a computation that has been constructed in a 
quasi automated way interfacing the GoSam and MadDipole programs. We calculate the 
electroweak NLO corrections to the pr oduc tion of a W + boson in association with two jets, 


which has also been calculated in Ref. 


138]. In the present calculation we treat the W-boson 


as a stable particle. Although this might be of limited interest from a purely phenomenological 
point of view, the decay of the W does not add any fundamental complication to the 
calculation which is not already present in the approach of a stable W-boson. Of course, 
in the case of an off-shell W production, the calculation should be rephrased in the context 
of the complex mass scheme; moreover, the IR subtraction procedure might be modified in 


order to have an efficient in tegra 


the presence of resonances 


215, 


don of real matrix elements plus local counterterms also in 


216]. The complete NLO computation should include the 


decay of the gauge boson together with the contributions of photon induced processes. For 
simplicity, we neglect both these contributions although they could be computed in the same 
framework of our calculation. By the way, these contributions will not affect significantly 
the inclusive results, but in particular the photon initiated processes can be important in the 
tails of distributions: indeed, they could well be the subject of a dedicated phenomenological 
study |244 ]. For all these reasons, we stress that in this review the process W + 2 jets serves 
as a proof of concept where one faces many subtleties of an electroweak computation while at 
the same time minimizing the necessary computational efforts. I 11 fact, in many aspects this 
process is identical to the calculation of the dijet process, which has been presented in Ref. 119 ], 


A. Computational setup 

The starting point is a Born contribution of the order 0(a^a). At this order in the coupling 
constants two classes of diagrams contribute to the Born process: the ones with two external 
quarks (with two additional gluons either radiated off the quarks or produced via gluon splitting) 
and the ones with four quarks involving a gluon exchange between the two fermionic currents. 
In principle there would be also a Born contribution of the order 0(a 3 ) involving only diagrams 
with four quarks where the internal gluon is replaced by a photon or a massive gauge boson. 
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FIG. 3: Example diagrams for a QCD Born (left diagram) and an electroweak Born (right diagram). 
At leading order only the left diagrams are taken into account. 

However, for the NLO EW corrections to W+2 jets, one can safely neglect the 0(a 3 ) underlying 
Born: it would lead to 0(a 4 ) contributions which will be much smaller than the 0(a 2 a 2 ) ones 
due to the different sizes of the coupling constants. Two sample diagrams are depicted in Fig. [3] 


The fact that the bh-boson is treated as a stable particle raises an additional complication 
at one loop. In fact, the building blocks of the calculation involve diagrams with internal 
massive propagators that can give rise to singularities when the momentum flowing in these 
propagators approaches the on-shell limit. On one hand these singularities could be avoided by 
means of the complex mass scheme, which would imply that the on-shell prescription for the 
process W + 2 jets should be abandoned. On the other hand, a naive inclusion of the widths 
in the massive propagators would lead to inconsistencies because of the different treatment of 
the external and the internal W bosons (in particular, the soft singularities corresponding to 
the W — 7 splitting would be turned into logarithms of the widths). However, in the process 
W + 2 jets at 0(a 2 a 2 ) these potentially problematic diagrams are interfered with a QCD Born 
that does not involve internal massive particles. This means that these terms do not give 
rise to physical resonances, however they are only integrable as principal value integrals. For 
a numerical Monte Carlo integration this is however problematic as it would simply try to 
integrate over a singula r reg ion. 


As discussed in Ref. 


138 ], one possible solution is to introduce a regulator width T in those 


pseudo-resonant diagrams that are not associated with IR singularities, so that the remaining 
diagrams are free of pseudo-resonances. Though this procedure technically breaks gauge invari¬ 
ance, the violation of gauge invariance is of order T /M. For the one loop contribution we adopt 
a different method that is based on a local K-factor. For each phase space point we require 
that the finite part of the one loop contribution normalized to the underlying QCD Born (i.e. 
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a Born that is free of resonances) has to be smaller than a certain value. The advantage of 
this approach is that one does not have to introduce a width for only specific diagrams. Nev¬ 
ertheless it can be related to the method of using a regulator width: in fact, a large K-factor 
is obtained if the internal particle is close to being on-shell and the K-factor effectively cuts 
the contribution off if the particle is too close to being on-shell. One can roughly estimate the 
width corresponding to a given K-factor by comparing the massive propagator to a massless 
one, requiring that 


p 


<K, 


(41) 


p* — m* 

and solving Eq. (PTTT) for p for a given value of K . This way one obtains the minimally allowed 
distance from the mass m. It is worth noting that the procedure in Eq. fTXTjl is only an approxi¬ 
mation, as in general many other aspects enter in the K-factor. As a default we used a value of 
1000 for the K-factor, but we varied its value from one hundred to one thousand and found that 
the results for the different values of the K-factor agree within their statistical errors. To further 
reduce numerical instabilities we generated the diagrams with a Higgs exchange separately and 
used an integration channel that resembles exactly the peak structure of these diagrams. This 
may serve as an additional strategy to improve numerical convergence for parton level results: 
loop diagrams from which one expects different peak structures can be generated separately 
and then integrated with a dedicated parametrisation of the phase space. 

For the subtraction terms consisting of a gluon exchange between an EW Born and a QCD 
Born as illustrated in Fig. [I] one can in principle use the same method, which is however not 
very convenient, as one has to perform this check for each subtraction term with its remapped 
kinematics. And it also requires to calc ulate the matrix element of the QCD amplitude. There¬ 


fore we adopt the method of Ref. 


138] and used a regulating width for the internal W- and 


Z-bosons. We varied this with between 0.1 and 1 GeV and found agreement between the results 
within their statistical uncertainties. 


1. Virtual corrections 

For this process we have to distinguish two different types of virtual corrections. Starting 
from the QCD like born structure we have to calculate the one-loop electroweak corrections 
leading to a contribution of 0(a 2 a 2 ). However this order in perturbation theory can also 
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be achieved by interfering a loop amplitude of 0[c? s a x ^ 2 ^ with a pure electroweak tree-level 
amplitude of 0(a 3 / 2 ). This contribution can be understood as an interference term between 
an electroweak tree-level amplitude with a QCD one-loop amplitude originating from the 
QCD born. The two different contributions are illustrated in Fig. [5j The lower two plots 
show an example diagram for the two cases. The left diagram contains a loop with both a 
gluon and a photon (or a weak gauge boson) in the loop. This is a good example that the 
distinction between QCD corrections and electroweak corrections does not make sense here, as 
this diagram can be regarded to be a QCD correction to an electroweak Born as well as an 
electroweak correction to a QCD born. This also means that these kinds of diagrams might 
contain singularities from both interactions, QCD and electroweak. The loop on the right 
hand is a pure QCD correction and therefore only contains QCD singularities. 


Simply speaking this situation always arises in the case where at tree-level the desired final 
state can be obtained in two different ways, either with a QCD Born or with an electroweak 
Born. The two contributions can be calculated independently, however, only the sum of the 
two yields a physically meaningful result and one has to sum up all the terms that contribute 
to a given order in perturbation theory at next-to-leading order, both in a s and a. 

There is one subtle difference between the two contributions. In the first case one starts 
from a QCD born amplitude with a given order in a and a s and the one-loop amplitude just 
adds another factor of a while keeping the original order in a s . This is also the situation 
one encounters performing QCD corrections and this is in some sense the standard way. The 
second contribution however is different. Starting from an electroweak born amplitude, the 
corresponding one-loop amplitude is two powers higher in a s but also one power lower in a. 
In view of automating electroweak calculations and combining different tools, in particular 
one-loop providers (OLP) and Monte Carlo tools (MC), this is a special case that is not 


included in the currently used standard interface 


217 


218]. Therefore an appropriate extension 


of the existing interface would be highly desirable. 


In this calculation the virtual amplitudes are generated with GoSam 


HEi. 


The GoSam 


framework is based on an algebraic generation of D -dimensional in 
diagrammatic approach, employing QGRAF [219 ] and FORM 


220 


egrands using a Feynman 


221 ] for the diagram gen- 
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FIG. 4: Schematical representation of the two different kinds of insertion operators, representing QED 
singularities (l.h.s.) and QCD singularities (r.h.s.). This also illustrates the two different types of real 
emission contributions. 


eration, and Spinney (2221 ]. Haggies [2231 ] and FORM 


;o w rite an optimized Fortran output. 


For the reduction of the tensor integrals we used Ninja 


224 


226], a package for the integrand 


reduction via Laurent expansion. One can however also use different reduction techniques such 


as integrand reduction via the OPP method 


227, 


228] as implemented in SAMURAI_15] 


or one can use methods of tensor reduction as the ones contained in Golem95 


The remaining scalar integrals have been evaluated using OneLoop 233 ]. All these methods 


229-232 . 


and tools have originally been developed for the case of QCD one-loop corrections. However 
for none of them it matters whether one uses them for QCD corrections or for electroweak 
corrections. The generation of the diagrams, and the reduction of the tensor integrals or the 
integrand reduction respectively rely on general properties that do not depend on the type of 
the interaction. 

As GoSam allows to specify the orders of tree-level amplitude and one-loop amplitude sepa¬ 
rately, i.e. does not rely on the fact that the one-loop amplitude is just one order higher in the 
coupling constant, it is straightforward to obtain the two different virtual contributions. We 
have extended the GoSam framework to incorporate the electroweak counter terms as well as 
the electroweak infrared insertion operator. 

The advantage of the inclusion of the insertion operator is twofold. In GoSam the agreement 
of the poles is the first level of the trigger for numerical instabilities. The numerical structure 
of the insertion operator is very simple and it is proportional to the Born. However on the side 
of the one-loop amplitude this simple results is obtained in a highly non-trivial way. Numerical 
instabilities can be detected as deviations from the correct coefficient of the pole terms. This can 
then either be used to trigger some kind of rescue system or it can be used to make quantitative 
estimation on the accuracy of the obtained result. As the pole cancellation only works for the 
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FIG. 5: Different virtual contributions. In the upper row the left plot denotes the interference between 
a QCD tree-level diagram with an electroweak one-loop diagram, the right plot displays the interference 
between an electroweak tree-level diagram and a QCD one-loop diagram. As the two pieces contribute 
at the same order both have to be taken into account. The lower row shows an example diagram for 
each of the two types of one-loop amplitudes. 

renormalized amplitude this is also an important check for the electroweak counter terms, 
which are new in the GoSam framework. Especially for these types of processes, where one 
has different virtual contributions, the check of this pole cancellation is very important as the 
cancellation only takes place when combining electroweak and QCD one loop contributions. 
We should stress that the formulas for the insertion operator in Eq. (12T1) are only valid for 
contributions from the l.li.s of Fig. [5j They are obtained by photon insertion between the 
underlying QCD amplitudes, i.e. the insertion operator only covers the QED singularities of 
these diagrams. It does not cover any singularity of QCD origin which are contained in the 
diagrams of the r.h.s. of Fig. 0 Therefore, in general one does not find pole cancellation between 
the two terms in Fig. 0and the QED insertion operator. The singularities from contributions 
stemming from the r.h.s of Fig. [5] are of purely QCD origin. This means that we need to 
apply a QCD insertion operator. In contrary to normal QCD corrections, the born-like matrix 
element is given by the interference term between a QCD Born amplitude and an electroweak 
Born amplitude with an additional gluon insertion that modifies the color structure of the 
interference term (as it does in normal QCD). The singularities from this insertion operator 
will then cancel against the QCD singularities from one-loop contributions of both l.h.s. and 
r.h.s diagrams in Fig. [5j The two different types of insertion operators are illustrated in Fig. |4j 
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For the production of a IF-boson in association with jets the number of processes that need to 
be generated can be reduced to: 

uu —> W + ud 
ud —> W + cc 
ud —> W + ss 
ud W + gg 
uu —> W + cs 
dd —> W + cs 

ud -> W + dd. (42) 

Note that for the process ud —» W + gg there exists no underlying electroweak born, hence this 
process does not have a QCD one-loop contribution. 

We have checked the pole cancellation mentioned above for all the independent subprocesses. 
All other processes can be obtained by crossing. Besides the renormalization of the electroweak 
sector this type of process also requires the renormalization of QCD singularities. They arise 
from the QCD loop diagrams that are interfered with an electroweak born process. The renor¬ 
malization procedure itself is straightforward and identical to a conventional QCD calculation, 
however the interpretation is a bit different. The underlying born that needs to be renormal¬ 
ized is an interference term between the QCD born and the electroweak born rather than the 
absolute square of an amplitude at a given order in a s . 

2. Regularization scheme dependence 

In the GoSam framework dimensional reduction is the more appealing approach as one deals 
with a simplified D-dimensional Dirac algebra and does not need to take care about additional 
renormalization terms for 75 . This typically leads to more compact expressions and hence is 
the method of choice. It is worth mentioning that for the given process, the scheme dependence 
for the QCD loops is proportional to the interference term between the QCD Born and the 
electroweak Born, with an additional modified color factor. For subprocesses that only have 
one color structure the scheme dependence vanishes. This is however not a general feature and 
does not hold for more than one color factor. 


36 


3. Real emission 

The real emission contributions contain processes based on an underlying QCD Born with 
an additional photon. We neglect additional radiation of massive electroweak gauge bosons. 
As discussed above, a second contribution stems from the radiation of an additional gluon from 
an underlying electroweak Born amplitude, interfered with a QCD Born amplitude. This con¬ 
tribution is necessary to cancel the QCD like singularities from the virtual corrections. All real 
emission contributions and corresponding subt racti on terms for both QCD and QED radiation 


have been generated with MadGraph 4 


3, 


235 ] and M ADPipole 44-1461], containing sub¬ 


traction terms based on the dipole formalism as developed in 


153 


196 ], The QED version of 


MadDipole has been extended to cover the photon radiation off the W-boson. The integra¬ 


tion over tree-level and real emission phase space has been performed using MadEvent 


236], 


The inclusion of the virtual contributions has been obtained using the BLHA interface, which 
has been adapted to be usable for electroweak corrections. 

4- Parameters 

As a consistent set of input parameters we use M\y , M z and GQ with 


M w = 80.398 GeV, M z = 91.1876 GeV, G^ = 1.1663787 • 10“ 5 GeV , (43) 


while a is a derived quantity: 


= 


y/2G u M* 


w 


7T 


Mb 

Ml 


(44) 


The Gn input scheme avoids large logarithms involving light fermion masses originating from 
the running of a from the Thomson limit to the electroweak scale 1 19U17l| . As these logarithms 
are effectively included in the definition of the coupling constant, they have to be subtr acted 
from electric charge renormalization counterterm with the following replacement 

1 


49 


119 


1711 ]: 


5Z e \ Gu = 5Z e 


l“(0) 


-Ar, 


(45) 


where 5Z e 


decay 


183 


□HU 


is defined in Eq. (HU while Ar summarizes the radiative corrections to muon 


185] apart from the QED corrections which coincide with those of the Fermi model. 
Due to the real emission of the photon, which happens at q 2 = 0 one power of the coupling 
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constant a is taken to be as a = a(0) = 1/137.035999074. Moreover, the one loop renormal¬ 
ization is a function of the Higgs and the top quark mass. In our numerical studies we used 
the values: 


m e = 510.99892 KeV, = 105.658369 MeV, m T = 1.77699 GeV, 


m t = 172.9 GeV, 


M h = 120 GeV. 


(46) 


As the W is treated as a stable particle in this calculation, no complex mass scheme is applied 
here. 

Renormalization and factorization scale are both set to Mw, the center of mass energy 
is chosen to be 14 TeV. As this is a calculation which is leading order in QCD we use the 
MSTW20081o pdf set. The CKM matrix has been chosen to be a diagonal unit matrix. In order 
to set cuts in an infrared safe way we use the following approach of democratic clustering: all 
particles, including the I F-boson a re passed through the jet algorithm and clustered according 


to an anti-kt algorithm 


233, 238] contained in the Fastjet package 239 ], In a second step 


we however require that the IF-boson as well as two jets containing a QCD parton are tagged. 
This ensures that we are also infrared safe from a QCD point of view. The resulting jets need 
to fulfill the following basic cuts: 


p T > 20 GeV, \r)\ < 4.4, A R > 0.4 . 


(47) 


This is in particular also true for the IF-boson. 


B. Numerical Results 

With the set of parameters stated above we find a total cross section at leading order of 

&tot,LO = 10873 ± 6pb . (48) 

As expected, the electroweak NLO corrections only have a small impact on the total cross 
section. Their contribution is negative and we find a total cross section of 

&tot,NLO = 10789 ± 6pb . (49) 

Therefore we find a correction of ~ 1% to the total cross section. The errors in both eqs. (1481) 
and (l49jl represent the 1 o Monte Carlo integration error. 
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FIG. 6: Differential distribution of the transverse momentum of the VF-boson (left plot) and the Ht 
distribution (right plot). 

More interesting than the total cross section are however differential distributions as one 
expe cts deviations stemming from contributions of large logarithms of the W mass (see e.g. 


240]). In Fig. |6] we show the transverse momentum of the W boson. Although the impact 


in the low energy region is negligible we do not obtain a constant K-factor but we see a 
substantial decrease of the cross section in the region where the transverse momentum of the 
W is large compared to its mass. In this region electroweak corrections typically supersede 
QCD corrections leading to a ~ 20% correction in the region of 500 GeV. A similar behavior 
is observed for the Ht distribution, also leading to a ~ 20% in the high energy tail of the 
distribution. This behavior stresses the importance of electroweak corrections at the LHC in 
the context of searches for new physics, as one would typically search for new physics in the 
high energy limit. For observables that are sensitive to Sudakov logarithms, like the transverse 
momenta or the Ht distribution, the size of the electroweak corrections grows beyond the 
uncertainty from the QCD corrections in the tails of the distributions. 

An interesting behavior can also be observed for the difference in rapidity between the W 
and the jets A yj\,w- Fig. [Tjshows this observable for the leading jet, i.e. the jet with the highest 
transverse momentum. One observes a negative correction for small differences in rapidity, but 
a positive correction for large values of A Uji t w which are however not relevant for the total cross 
section, as the contribution to the total cross section from this region is small. The radiation of 
a particle with a high transverse momentum will more likely happen to take place in the central 
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FIG. 7: Differential distribution for absolute value of the rapidity difference between the W and the 
leading jet. 

region, as one needs less energy to get the same transverse momentum as in a more longitudinal 
direction. A bF-boson with a high transverse momentum has to recoil against the jets (or at 
least the leading jet). Therefore we expect the events with a high transverse momentum of 
the W happening for small differences in rapidity, which explains the negative next-to-leading 
order corrections in this region. Turning this argument around we expect events with large 
differences in rapidity to come with smaller values of the transverse momentum. This leaves 
more phase space to additional radiation which comes with a positive contribution and hence 
shifting the NLO corrections to positive values. 

The same behavior is also found for the rapidity distribution of the W which is shown 
in the left plot of Fig. [SJ although the rapidity distribution to a larger extents mixes events 
with different transverse momenta and therefore the effect is much less pronounced. Generally 
speaking we do not expect visible effects of electroweak NLO corrections for observables where 
the events with high transverse momentum of the W are evenly distributed as a function of the 
observable under consideration. One such example is shown in the right plot of Fig. [8l where 
the difference in the azimuthal angle between the two leading jets is displayed. In that case we 
obtain a flat K-factor yielding a small negative correction at next-to-leading order. 








40 




FIG. 8: Differential distribution of the rapidity of the VF-boson (left plot) and the difference of the 
azimuthal angles of the two leading jets (right plot). 

C. Virtual amplitudes for single phase space point 


In this section we give the results for all the ingredients that enter the virtual contribution 
for a single phase space point reported in Table I. The regularization scheme is DRED, the 
couplings g s and e are set to one for simplicity and all the results in the tables do not contain 
the overall factor a^ s )/2n. The other parameters are set as stated above. The notation is as 
follows: Born-QCD denotes the tree-level matrix element squared that is of the order 0(a^a), 
while Born-EW denotes the tree-level matrix element squared of the order 0(a 3 ) |245|. The 
first line in the tables gives the contribution of the unrenormalized clectroweak corrections, 
the second line the corresponding counter terms. The third line denotes the unrenormalized 
contribution of QCD-loops interfered with the electroweak tree-level amplitudes, followed 
by QCD renormalization terms. The last two lines then denote the contributions from the 
electroweak and QCD insertion operators. 
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E 

Px 

Py 

Pz 

250.0 

0.0 

0.0 

250.0 

250.0 

0.0 

0.0 

-250.0 


111.48897463938037 56.330449672673069 -24.583961957306506 46.781003797444882 
207.38546864044983 -167.21520572166636 -118.01342209209260 -33.475958310305167 
181.12555672016978 110.88475604899321 142.59738404939910 -13.305045487139747 


TABLE I: Phase space point where all the relevant amplitudes of Eq. (1421) have been evaluated. 

1. uu —>• W + ud 

Born-QCD: 0.3331176339143377E - 03 (Born-EW: 0.1159155705277286E - 02) 


Process 


uu —»■ W + ud 



l/e 2 

l/e 

finite 

EW-Loop unren. 

0.4515129528621876E-03 

0.5179344797596786E-02 

-0.4048755402433699E-02 

EW-Ren. 

0 

-0.4682725341929206E-02 

-0.6484221636837220E-03 

QCD-Loop unren. 

0.9326828685162938E-03 

-0.2347576391436467E-02 

0.1872325091832032E-02 

QCD-Ren. 

0 

0.1224146264927649E-02 

0.8851406118179909E-05 

I-EW 

0.4811699156540433E-03 

-0.1471889938882662E-02 

- 

I-QCD 

-0.1865365737032609E-02 

0.2098700609724008E-02 

- 

E 

-0.8E-16 

0.1E-15 

-0.2815845104231514E-02 
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2. ud —>• W + cc 


Born-QCD: 0.1509105387580410E - 03 (Born-EW:0.3226028869165598E - 03) 


Process 


ud —>• W + cc 



l/e 2 

l/e 

finite 

EW-Loop unren. 

-0.2179818893171668E-03 

0.2371604818917693E-02 

-0.7006797514364052E-03 

EW-Ren. 

0 

-0.2187116288266569E-02 

-0.3231845769861692E-03 

QCD-Loop unren. 

0 

-0.7013028153765194E-05 

0.1644192724679622E-04 

QCD-Ren. 

0 

0 

0 

I-EW 

0.2179818893171704E-03 

-0.1915015588048693E-03 

- 

I-QCD 

0 

0.14026056307548736e-04 

- 

E 

0.4E-17 

0.4E-16 

-0.1007351745650022E-02 


3. ud —>• W + ss 


Born-QCD: 0.1509105387580410E - 03 (Born-EW: 0.1683155172046103E - 03) 


Process 


ud —> W + ss 



l/e 2 

l/e 

finite 

EW-Loop unren. 

-0.1173748634784783E-03 

0.2359759858112650E-02 

-0.8847947525100161E-03 

EW-Ren. 

0 

-0.2201404574690460E-02 

-0.3267300831621827E-03 

QCD-Loop unren. 

0 

0.1073376195468900E-04 

-0.8486634366399400E-04 

QCD-Ren. 

0 

0 

0 

I-EW 

0.1173748634784764E-03 

-0.1476215214675483E-03 

- 

I-QCD 

0 

-0.2146752390936322e-04 

- 

E 

-0.2E-17 

-0.3E-16 

-0.1296320523810437E-02 
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4- ud —>• W + gg 


Born-QCD: 0.4421271590046107E - 03 (Born-EW: 0) 


Process 


ud —> W + gg 



l/e 2 

l/e 

finite 

EW-Loop unren. 

-0.2456261994470038E-03 

0.5954936826411168E-02 

-0.4059861036705788E-02 

EW-Ren. 

0 

-0.5796026153495491E-02 

-0.6841797425946791E-03 

QCD-Loop unren. 

0 

0 

0 

QCD-Ren. 

0 

0 

0 

I-EW 

0.2456261994470059E-03 

-0.1589106729178175E-03 

- 

I-QCD 

0 

0 

- 

E 

0.2E-17 

-0.2E-14 

-0.4218771709623606E-02 


5. uu —>■ W + cs 


Born-QCD: 0.5471920731180749E - 04 (Born-EW: 0.1560261025155393E - 03) 


Process 


uu —» W + cs 



l/e 2 

l/e 

finite 

EW-Loop unren. 

-0.7903885500595288E-04 

0.1039825482424979E-02 

0.6729212743821051E-03 

EW-Ren. 

0 

-0.7933882016361721E-03 

-0.1173430614024200E-03 

QCD-Loop unren. 

0 

0.4658874285701936E-05 

-0.6180475632648689E-04 

QCD-Ren. 

0 

0 

0 

I-EW 

0.7903885500594414E-04 

-0.2417784065030774E-03 

- 

I-QCD 

0 

-0.9317748571408931e-05 

- 

E 

-0.9E-17 

0.2E-16 

0.49379907590015959E-03 














44 


6. dd —>• W + cs 


Born-QCD: 0.5471920731180750E - 04 (Born-EW: 0.8665998896761033E - 03) 


Process 


dd -> W+cs 



l/e 2 

l/e 

finite 

EW-Loop unren. 

-0.4255938346473633E-04 

0.8312276950839127E-03 

0.1195260305972849E-02 

EW-Ren. 

0 

-0.7986225826314173E-03 

-0.1186419246527548E-03 

QCD-Loop unren. 

0 

-0.3536873862079221E-05 

-0.2058866127333486E-03 

QCD-Ren. 

0 

0 

0 

I-EW 

0.42559383464739I6E-04 

-0.3614198631460137E-04 

- 

I-QCD 

0 

0.7073747724144110E-05 

- 

E 

0.2E-17 

-0.4E-16 

0.8707573878337068E-03 


7. ud -> W+dd 


Born-QCD: 0.8493852813137004E - 03 (Born-EW: 0.1939842768622516E - 02) 


Process 


ud -> W + dd 



l/e 2 

l/e 

finite 

EW-Loop unren. 

0.1913723183724379E-02 

0.9838014169469896E-02 

-0.1220570133146914E-01 

EW-Ren. 

0 

-0.1197778947261055E-01 

-0.1665029533330621E-02 

QCD-Loop unren. 

0.2574356180301684E-02 

-0.6257504923291653E-02 

0.1973036296675117E-02 

QCD-Ren. 

0 

0.3378842486646013E-02 

0.2443131831181329E-04 

I-EW 

0.6606329965773225E-03 

-0.8308733675698238E-03 

- 

I-QCD 

-0.5148712360603448E-02 

0.5849311107356306E-02 

- 

E 

-0.6E-16 

0.2E-15 

-0.1187286557206206E-01 


V. CONCLUSIONS 


The recent start-up of the LHC’s Run 2 will allow to test the Standard Model through a 
set of new precision measurements. The search for new physics beyond the Standard Model 
is just started and new regions of phase space will be explored. The effect of the new physics 
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might be large, but given the unique opportunity offered by the LHC physics program one has 
to look also for small deviations from the Standard Model. From the theoretical point of view 
this requires predictions based on precise and reliable computations whenever possible. This 
triggered the so called NLO revolution that was originally focused on providing NLO QCD 
corrections in an automated way. For many important processes like Higgs and top-quark 
production for instance the NLO corrections were not precise enough to meet the requirements. 
This lead to an ongoing active field of research namely pushing the limit further to NNLO 
QCD corrections. On the other hand also electroweak corrections will play an important role 
in a large number of searches. Although their impact on the total cross sections is typically 
rather small, they can lead to significant differences in distributions, in particular in the high 
energy tails. However it is exactly the high energy limit that is sensitive to new physics effects. 
Therefore, for a reliable comparison with LHC data from Run 2, the inclusion of electroweak 
higher order effects becomes inevitable for a large class of processes (also with large final state 
multiplicities). In this context, it is clear that fully automated tools for the computation of 
NLO electroweak corrections will play a crucial role at the LHC and at future hadron colliders. 

We have presented an overview on the current status and recent development in the held 
of electroweak NLO corrections for LHC physics. In particular we focused on the topic of the 
automation of O(a) calculations, which is a very active and fast developing area of research. 
We have summarized the existing tools and the phenomenological predictions that have been 
obtained for LHC relevant processes. Furthermore we have discussed some of the general issues, 
related to the calculation of either the virtual one loop corrections or the real emission contri¬ 
butions, that have to be faced in order to develop an automated tool for the computation of 
NLO electroweak corrections starting from already available programs such as GoSam (for the 
generation of one loop amplitudes) and MadDipole (for the subtraction of infrared singular¬ 
ities). We discussed the regularization scheme dependence of the virtual one loop amplitudes, 
the renormalization counterterms and infrared subtraction terms (at least for the specific case 
of the CDR and DRED schemes) and we provided the corresponding scheme transition rules 
that are required in order to add up consistently these three building blocks of the calculation. 
We also pointed out the 11011 trivial aspects of the infrared subtraction procedure in presence of 
QCD and QED-like infrared singularities. Lastly, as a non-trivial example, we have presented 
the full electroweak NLO corrections to the production of a W-boson in association with two 
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jets within the GoSam+MadDipole framework. 
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